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The random percolation model can be viewed as the dual of a well defined confining gauge theory; 
since this theory, having no Monte Carlo dynamics at all, is simple to simulate, it is possible to 
study the properties of the flux tube with very high precision; we show it can be described by 
the effective string picture. Our results are lattice regularisation independent, therefore they are 
well defined also in the continuum limit, and, for the first time in a gauge theory, it has been 
possible to determine the next-to-leading quantum corrections throughout the computation of the 
T^ coefficient of the Taylor expansion of a{T). Furthermore, this coefficient results to be related 
to the universal ratio T^ y/oQ. 
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1. Introduction: the effective string theory 

The assumption behind the effective string theory is that the color flux connecting a pair of 
quark is squeezed inside a thin flux tube; as a consequence, the confining potential rises linearly. 
According to this picture the flux tube should behave as a free vibrating string. 

Unfortunately, the action of this effective theory is unknown; the simplest assumpion is that it 
is described by the Nambu-Goto action, which is proportional to the string worldsheet area. 

We try to summarize the outcome of many studies on this argument discussing some properties 
of the first terms of the low temperature expansion of the string tension (the first term Oq is the zero- 
temperature string tension): 

a(r) = ao-(^/-2)^r2+£c„r . (1.1) 

The second term, the analogue of the Liischer term at finite temperature which was calculated in 
Ref. [[T]], does not depend on the gauge group and is expected to be independent of the interaction 
terms of the effective theory. Thanks to a certain open-closed string duality it was shown that for 
any dimensionality C3 = and, in three dimensions, C4 is universal ; hence, it coincides with the 
value calculated in the NG model [Q, ^: 

C4 = -(^/-2)-^. (1.2) 
72 ao 

Using a different approach to the effective string theory, Ref. [Q], the above results were confirmed 
for all values of d. 

In this paper we will evaluate the coefficients c„ up to n = 6, in a simple, but not trivial, model: 
the gauge theory dual to the 3d percolation model. All of the date agree with the universal values 
of C2 and C4 and lead to C5 = and c^ = %^ j (Cag) , where C ~ 300. 

We decided to focus our attention to the behaviour of the Polyakov-Polyakov correlation func- 
tion {P(0)P*{r)) at finite temperature 7 = ^ in (2+ 1) -dimensions; r is the distance between the 
Polyakov loops, £ is the time extent of the lattice and a is the lattice spacing. 

The functional form of the correlator has been calculated at the next-to-leading order (NLO) 
in Ref. |@1: 

(P(0)P*(r)) _ e- {d-2)nH[2E,{T)-Ei{T)] . ^ ( I \\ 
{P{0)P {r)Wo - (I Yj^, + ^ j j ' ^ 

where 77 is the Dedekind function, E4 and E2 the two Eisenstein functions and z = ii/ (2r). 



Using Eq. (1.3) one can find, for asymptotically large r: 



a(r) = a - - = ao - fr^ - + 0{T') . (1.4) 



6 72a 6 72ao 

Note here the difference between a and ao: a = ao + 0{T^) 
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2. The gauge theory dual to the percolation model 

In this paper we study a particular gauge theory, first introduced in Ref. [^], that is dual to the 
random percolation model. A more complete account will be presented in Ref. [^]. 

It is known in three dimensions it is possible to study a well defined Sg-gauge' theory dual 
to the 2-state Potts model through the Kramers-Wannier duality [Q]. Thanks to the fact one can 
map some gauge invariant observables, such as Wilson loops and Polyakov correlators, into the 
corresponding quantities of the spin model, it is numerically convenient to inspect the properties of 
the dual theory instead of those of the gauge model. The ingredient which is then fundamental in 
our approach is the Fortuin-Kasteleyn reformulation |^] of the Q-state Potts model, by which it is 
possible to determine gauge observables in a very efficient way. This approach, that for 2 > 1 is 
only a powerful numerical method, can be applied to the random percolation model whose gauge 
formulation is not known: it is the gauge theory in the I limit. 

The key ingredient is the method used to calculate the Wilson loops in this setup: we define a 
procedure to determine its value studying some topological properties of the dual model. 

The connected components of the graph, formed by active links, are known as clusters; Wy is 
the value of the Wilson loop associated with a loop with contour y. We set = 1 if there is no 
cluster topologically linked to the contour y, otherwise we set Wj = 0. The same linking properties 
are used to determine the Polyakov-Polyakov loop correlators {P{0)P*{r)): at finite temperature 
the contour 7 is a r x £ rectangle, with two sides identified. 

Another interesting study of this model, related to the monopole mass, can be found in [|l^]. 

3. Simulations 



The idea behind this work is not only to verify whether one can observe the presence of shape 
effect due to rough fluctuations of the string, in agreement with the universality predictions of the 



effective string picture (as a matter of fact we have discussed this point in Ref. Jllp); we also would 
verify that our results are not regularisation dependent. In other words, we would discuss if our 
results describe a "real" phenomenon and not a lattice/model artifact. 

We therefore study five different systems (see Ta- 
ble [1]) characterized by different occupancy probabil- 
ity p, different kind of percolation (site or bond) and 
different geometry of the lattice (simple cubic lattice 
(SC) and body-centered cubic lattice (BCC)). 

We worked on a lattice of size L? x £, where i, 
the inverse of the temperature, was chosen such that 
0-371 < r < O.STc. The value of the spatial size was 
L = 128 which was in most cases sufficient to account 

for the infinite volume limit. Just in the case ic = 8, simulated at £ = 10 and £ = 11, we found a 
sizable dependence on the lattice size L; in this case we performed further simulations on larger 
lattices in order to extract the correct value of a using the scaling relation (v = 4/3 is the termal 
exponent of 2d percolation model): 

-cL-i/V (3.1) 



Lattice 


P 


4 = l/aTc 


SC bond 


0.272380 


6 


SC bond 


0.268459 


1 


SC bond 


0.265615 


8 


SC site 


0.3459514 


7 


BCC bond 


0.21113018 


3 



Table 1: The five systems simulated. 



Sq is the symmetric group. 



3 



Confining string beyond the free approximation: the case of random percolation 



Pietro Giudice 



For each system, we measured ( P{0)P* (r) ) by varying the distance between the two Polyakov 
lines from r = 8 to r = 50; to reach an acceptable statistics, we collected data from 10^ configura- 
tions for each value of p and £. 

The algorithm used, described in detail in Ref. [^, is basically aimed at determining the linking 
properties of clusters with the Polyakov-Polyakov contour. 



4. Numerical results 



Our numerical results are compared with the expected behaviour of the Polyakov-Polyakov 
correlation function given in Eq. Being an expression valid in the infrared limit we use a 

sliding window analysis to determine the correct values of the fitted parameter a: we fitted the data 
in the range r,„,„ < r < r^^ix by progressively discarding the short distance data, varying rmin but 
fixing the value of Vmax = 50a: (see Fig. |]). In all five sistems considered a large plateau appears for 
all values of i not too close to Ic, showing the stability of the fits and so the suitability of the string 
picture to describe our data. It is important to note that, as Fig. |l| shows, there are different values 
of the string tension for different values of i, i.e. of T . In other words, the value of a is not yet the 



string tension at zero-temperature Oq, the formula Eq. ( |1.3D is not the exact formula because it only 
takes into account the temperature dependence up to the order (see Eq. (L4)). We studied the 
dependence of a on ^ and we verified, in all cases, that fox aT = \/i low enough the correction is 
proportional to (see Fig. Therefore, we used the value of a to determine the value of o{T) 



by Eq. (1.4), i.e. we reconstructed the correct dependence of the string tension on the temperature; 
then we used these data to perform a new fit to determine the first model-dependent term by means 
of the Ansatz: 

a{T) = ao-^T^-SzT^ + ^T'' + 0{T^) . (4.1) 



72a 



Col 



Thereby, we can identify stable values both for the zero-temperature string tension Oq and the 
coefficient C, see Table ^. 



Lattice 


4 = l/aTc 


C 






SC bond 


6 


291(7) 


0.012612(6) 


1.4841(4) 


SC bond 


7 


281(5) 


0.009234(5) 


1.4866(5) 


SC bond 


8 


297(5) 


0.007059(5) 


1.4878(5) 


SC site 


7 


307(9) 


0.009399(8) 


1.4735(6) 


BCC bond 


3 


295(14) 


0.0474(4) 


1.531(7) 



Table 2: The parameters C and a^ffo in the fit ( 4. 1 ) for the numerical experiments listed in Table |l]. The last 
column is the universal ratio Td as obtained by combining the second and the fourth columns. 

Note that the five values of C coincide up to the statistical errors. The value of / y/Oo, 
obtained by combining the precise determination of o^Oq with the deconfined temperaure Tc, is 
an important universal quantity which characterizes the particular gauge theory; the small varia- 
tions appearing in Table ^ are presumably due to the corrections-to-scaling that we have neglected. 
Nonetheless, we can assert the value closer to the continuum limit is that obtained in the simulation 
with bond percolation and l,. = 8 where statistical and systematic error were better under control; 
therefore we will use, in the following, Tcj ^/Oq = 1.4878(5). 
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If we plot the adimensional ratio o{T)/T} versus the reduced temperature t = {T — Tc)/Tc 
it turns out that all data lie almost on a unique curve, see Figure ^; this non-universal behaviour 
is related to the fact the five different systems are characterized by different "universal" value of 
the quantity T^/ ^/Gq and the adimensional variables used are very sensitive to it. As a matter of 
fact, if we impose the value of o{T)/T^ is the same for all systems, i.e. we determine for each 
system a new Tc value by which Tc/^/oq= 1.4878(5), all data fall on a unique universal curve as 
Figure ^ shows. This is the most important result of this work because it shows our results are 
independent of the regularisation used, therefore we are studing a "real" gauge theory well defined 
in the continuum limit. 

It is interesting to note that the values of C and of Tc/ ^/Oq can be determined with only two 
pieces of information: (1) the data are all in the scaling region and (2) they show a linear behaviour 
in the range —0.55 < t < —0.225. This means we can impose the two following equations to 
coincide in the above range (S = || and x = y)' 

_ c ^ 2 4 , 6 



S--x^ x^ + — ^.x"" , (4.2) 

6 725 CS^ ' 

Aix-l). (4.3) 



r^2 6 725 CS^ 



Immediately, without using numerical data, it is possible to determine C ~ 290 and Td — 
1.4884; these two values are remarkably close to those obtained using the numerical data, see 
Table ^ In Figure ^ we plot Eq. ( |4.2| ) (dashed line) and Eq. (13) (dotted line) using those values; 



the numerical data lie on the two curves in the scaling region. This is an important observation 
because it means the two quantity, C and r^/y^ are constrained to each other in that region. 



5. Conclusions 



In this paper we have studied, by numerical simulation, the gauge theory dual to the percolation 
model; we can conclude it is possible to describe the long distance dynamics of this theory by 
means of an effective string picture. Our numerical experiment demonstrate that the quantities 
which characterize the effective string theory do not depend on the specific regularisation used. 
Moreover, we determined with high precision the value of o{T)/T^ and, for the first time in a 
gauge theory, we have determined the value of the coefficient C of the string tension cj(r). 
Furthermore, it was possible to show that the universal ratio Tc/^/oq and the coefficient C are 
bound together in the scaling region. 
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Figure 1: The fitted value of the string tension a as a function of the minimal distance r,„,„ of the set of 
Polyakov-Polyakov correlators considered in the fit; case of bond percolation with 4 =8. 
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Figure 3: Plot of the scaling variable a{T)/T^ versus the reduced temperature. 
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Figure 4: Plot of the scaling variable a{T)/T^ versus the reduced temperature when we impose the value 
of Tc/y/a^ = 1.4878. Dashed hne is the plot of Eq. (4.2), dotted hne is that of Eq. (4.3). 



7 



